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A novel method to determine the density and temperature of a system is proposed based on 
quantum fluctuations typical of Fermions in the limit where the reached temperature T is small 
compared to the Fermi energy e/ at a given density p. Quadrupole and particle multiplicity fluc- 
tuations relations are derived in terms of — . This method is valid for infinite and finite fermionic 
systems, in particular we apply it to heavy ion collisions using the Constrained Molecular Dynam- 
ics (CoMD) approach which includes the Fermi statistics. A preliminary comparison to available 
experimental data is discussed as well. We stress the differences with methods based on classical 
approximations. The derived 'quantum' temperatures are systematically lower than the correspond- 
ing 'classical' ones. With the proposed method we may get important informations on the Equation 
of State (EOS) of quantum Fermi systems to order O(^j) 3 , in particular near the Liquid-Gas (LG) 
phase transition and at very low densities where quantum effects are dominant. 

PACS numbers: 42.50.Lc, 64.70.Tg, 25.70.Pq 



In recent years, the availability of heavy-ion ac- 
celerators which provide colliding nuclei from a few 
MeV/nucleon to GeV/nucleon and new and performing 
47T detectors, has fueled a field of research loosely referred 
to as Nuclear Fragmentation. The characteristics of the 
fragments produced depend on the beam energy and 
the target-projectile combinations which can be exter- 
nally controlled Q]-[|j|. Fragmentation experiments could 
provide informations about the nuclear matter proper- 
ties and constrain the EOS of nuclear matterQ. From 
conventional nuclear physics we know that there is a 
stable equilibrium state at the normal nuclear density 
pa = 0.145 — 0.17/m -3 with a compressibility in the 
range of K — 180 — 240 MeV and a binding energy of 
15-16 MeV/nucleon [H-[|. Even though a large variety of 
experimental data and refined microscopic models exist, 
to date it does not exist a method to determine densities 
and temperatures reached during the collisions, which 
takes into account the genuine quantum nature of the 
system. In this work we discuss some properties at fi- 
nite temperatures assuming either a classical gas or a 
quantum Fermi system. We show that at the densities 
and temperatures of interest the classical approximation 
is not valid. This is at variance with many experimen- 
tal and theoretical results in heavy ion collisions near 
the Fermi energy 0, @-[ll| which assume the classical ap- 
proximation to be valid. We base our method on fluctua- 
tions estimated from an event by event determination of 
fragments arising after the energetic collision. A similar 
method has recently been applied to observe suppression 
of fluctuations in a trapped Fermi gas[ll|- We go beyond 
the method of (l2| by including quadrupole fluctuations 
as well to have a direct measurement of densities and 
temperatures for subatomic systems for which it is dif- 
ficult to obtain such informations in a direct way. We 
also suggest a method for calculating an excitation en- 
ergy which should minimize collective effects and could 
be applied when a limited information is available, for 
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Figure 1. Temperature versus thermal energy per particle 
derived from quantum fluctuations (full symbols joined by 
dashed lines) compared to the classical case (open symbols). 
(Top) Circles refer to protons, squares to neutrons and tri- 
angles to protons and neutrons. (Bottom) Same as above 
for protons. Data: down triangles from classical quadrupole 
fluctuations [Sj], star symbols from particle ratios Q. 

example if only light cluster are measured. We apply the 
proposed method to microscopic CoMD approach [13| 
which includes Fermionic statistics. The resulting densi- 
ties and temperatures calculated using protons and neu- 
trons, even though surprising at first, give a perfectly 
reasonable EOS and a clear first order phase transition. 

A method for measuring the temperature was proposed 
in [9j based on momentum fluctuations of detected par- 
ticles. A quadrupole Q xy =< p x ~pi > is defined in a 
direction transverse to the beam axis (z-axis) and the av- 
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erage is performed, for a given particle type, over events. 
Such a quantity is zero in the center of mass of the equi- 
librated emitting source. Its variance is given by the 
simple formula: 



■i-y 



= d 3 P (pi-pirm 



(i) 



where f(p) is the momentum distribution of particles. In 
Q a classical Maxwell-Boltzmann distribution of parti- 
cles at temperature T c i was assumed which gives: a 2 y = 
NArr^T^, m is the mass of the fragment. N is the av- 
erage number of particles which could be conveniently 
normalized to one. In heavy ion collisions, the produced 
particles do not follow classical statistics thus the cor- 
rect distribution function must be used in eq.(l). Pro- 
tons(p), neutrons(n), tritium etc. follow the Fermi statis- 
tics while, deuterium, alpha etc., even though they are 
constituted of nucleons, should follow the Bose statistics. 
In this work we will concentrate on fermions only and in 
particular p and n which are abundantly produced in the 
collisions thus carrying important informations on the 
densities and temperatures reached. 

Using a Fermi-Dirac distribution f(p) and expanding 
to 0(T/e f )\ where e f = e/o(£) 2/3 = 36(^f/ 3 MeV is 
the Fermi energy of nuclear matter, we get |14j : 



xy 
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This result is in evident contrast with the classical one: 
even at zero T and ground state density po, quadrupole 
fluctuations arise from the Fermi motion, but those fluc- 
tuations cannot be observed since in this case the number 
of emitted particles N — 0. The quadrupole fluctuations 
depend on temperature and density through e / , thus we 
need more informations in order to be able to determine 
both quantities. 

Within the same framework we can calculate the fluc- 
tuations of the p,n multiplicity distributions. These are 
given by [l4j : 



< (AN) 2 > 
N 



2e/ £/ 



(3) 



This is a very simple result valid in the indicated approx- 
imation. The difference with the classical case is again 
striking (the ratio in eq.(3) equal to one for a classical 
perfect gas). This relation has also been derived and ap- 
plied to trapped Fermi gases in (l2^ . These quantities, 
eqs.(2-3), can be easily verified experimentally and the 
corresponding densities and temperatures can be evalu- 
ated for each physical situation. 

To illustrate the strength of our approach we simulated 
40 C a + 40 Ca heavy ion collisions at fixed impact param- 
eter b = lfm and beam energies Ei^/A ranging from 4 
MeV/A up to 100 MeV/A. Collisions were followed up 
to a maximum time t = 1000/m/c in order to accumu- 
late enough statistics. Particles emitted at later times 
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Figure 2. Classical temperatures vs quantum temperatures. 
Symbols as in figure 1, open symbols refer to Bauer's approx- 
imation, eq.(6). 



(evaporation) could affect somehow the results and this 
might be important especially at the lowest beam ener- 
gies. The choice of central collisions was dictated by the 
desire to obtain full equilibration. This however, did not 
occur especially at the highest beam energies due to a 
partial transparency for some events. For this reason the 
quadrupole in the transverse direction, eq.(l), was cho- 
sen. Furthermore, in order to correct for collective effects 
as much as possible, we defined a 'thermal' energy as: 



/Eth\ _ E cm -&p( n ) 3 .E p ( 



p(n)xy , 



p(n) 



(4) 



'-) are the average total and 



where (^) and 
transverse kinetic energies (per particle) of protons 



(and/or neutrons). Q va i 



N p 
Z{N) 



8MeV, similarly for 



protons plus neutrons. 8 MeV is the average binding en- 
ergy of a nucleon, Z (N) the total charge (neutron num- 
ber) of the system and N p ^ n j the average number of pro- 
tons (neutrons) emitted at each beam energy. Clearly for 
a completely equilibrated system all the center of mass 
energy, -^ys is converted into thermal energy (plus the 
Q value)- If n ot, our approximation will account for some 
corrections, and this will become more and more exact 
when many fragment types are included in eq.(4). How- 
ever, this approximation might be important in experi- 
ments where only some fragment types are detected or 
if, because of the time evolution of the system, different 
particles are sensitive to different excitation energies, for 
instance if some particles are produced early or late in 
the collision. 

In figure 1 (top) we plot the estimated temperatures 
at various 'thermal' energies both for the quantum (full 
symbols) and classical approximations (open symbols). 
As we see the quantum case is systematically lower than 
the classical one. We also notice a difference if the T are 
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estimated from the proton distributions (circles), or neu- 
trons (squares) or the sum of the two (triangles) . This is 
clearly a Coulomb effect which gets smaller as expected 
at higher energies as we will demonstrate more in detail 
below. The back-bending observed at T re 3MeV for all 
cases indicates a liquid-gas phase transition, in particular 
we observe that such a back-bending is more marked for 
the protons case as first discussed by Gross |15| . In the 
bottom part of figure 1, we compare the protons results to 
experimental data. The down triangles are derived using 
the 'classical' quadrupole fluctuations [9] thus should be 
very similar to our classical results and the agreement is 
reasonable at the lowest excitation energies. However, we 
stress that the experimental data were obtained for dif- 
ferent systems at a fixed 35MeV/A beam energy. In par- 
ticular projectile like fragments (PLF) were isolated and 
analyzed and the excitation energy was obtained from 
all fragments differently from eq.(4). Thus there might 
be a mismatch in the abscissa and this could be espe- 
cially important for large excitations. Also the detector 
acceptance might be important. Similar considerations 
apply to the data [8|| obtained using double particle ra- 
tios (star symbols) [fj. Even though in the latter case it 
is perfectly correct to use classical approximations since 
different particle types are used (usually a mixture of 
bosons and fermions @), the underlying assumption is 
that all those particles are sensitive to the same density 
and temperature. If T and p 'seen' by different particles 
are different, then the results give some kind of 'averag- 
ing' which might hide interesting quantum effects. In the 
top part of figure 1 we see that temperatures are differ- 
ent for protons and neutrons at a given excitation energy, 
thus we expect that other particles might give different T. 
This implies that different particle ratios might produce 
different results Q. 

Using eqs.(l-2), we can easily show that, in the region 
of validity, the 'classical' T c i is always larger than the 
'quantum' temperature T: 




(5) 



A similar result has been found by Bauer [l6| in order 
to explain the large 'a ppa rent' temperature observed in 
particles spectra. In [16j a relation between the final 
(classical) temperature T' cl and the input Fermi-Dirac T 
was found: 



5 12 e f 



(6) 



The ratio — entering the equations above can be directly 

obtained from eq.(3). Even though eqs.(5-6) might look 
different at first sight, they give very similar results as 
can be seen in figure 2 where the classical T c ; is plotted 
vs. the quantum one. Bauer's approximation, eq.(6), is 
given by the open symbols. Thus we argue that quan- 
tum temperatures are smaller than derived when fitting 
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Figure 3. Temperature divided the Fermi energy versus den- 
sity normalized to the ground state one derived from quantum 
fluctuations, eqs.(l-3). Symbols as in figure f. 



experimental results with a classical approximation. The 
reason of such small quantum temperatures is the Fermi 
energy entering eq.(5) or (6). In passing, we notice that 
our approximation, eq.(5), reproduces better the numer- 
ical results plotted in fig. 2 of ref . (To| . 

In figure 3 we plot the ratio 3L directly obtained from 
eq.(3), versus reduced density which is obtained from 
eqs.(2) and (3). The lowest T (highest j-) corresponds to 
the lowest beam energy as well and gives the lowest den- 
sity, especially for the neutrons case. This result might be 
surprising at first, but it simply tells us that at the lowest 
energies nucleons from the surface of the colliding nuclei 
come into contact. Those nucleons are located in a low 
density region, especially neutrons which do not feel the 
Coulomb field. With increasing beam energy, the over- 
lapping region increases and more and more fermions are 
emitted. At about 20 MeV/A a large number of nucle- 
ons are excited and the emission from surface becomes a 
volume emission. This explain the minimum in the plot, 
which is due to the increase of T and e/ when deeper 
regions of the nuclei are affected. Fragmentation starts 
around the beam energy which gives the minimum in the 
plot, where we observe a power law in the mass distribu- 
tion as well. The lowest density (as well as T) is explored 
by the neutrons only. It is important to stress that the 
ratio plotted in figure 3 is always smaller than one which 
confirms the approximations used in eqs.(l-6). 

The best way to visualize the results is by plotting the 
energy density e = {—fi L )p versus temperature as in figure 
4. Now different particle types scale especially at high 
T where Coulomb effects are expected to be small. A 
rapid variation of the energy density is observed around 
T re 2MeV for neutrons and T re 3MeV for the other 
cases which indicates a first order phase transitionfijj . 
Notice that a 'plateau' in the caloric curve i.e. (—V*-) 
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Figure 4. Energy density versus temperature. Symbols as in 
figure 1. 

vs T @, § has been experimentally observed around 
T c i s» 6MeV. Such a value agrees with our classical ap- 
proximation plotted in figure 1, but differs greatly with 
the quantum results, figure 4. It would be very interest- 
ing to reanalyze the experimental data in the framework 
of this paper. We also notice that Coulomb effects be- 
come negligible at around T — 3MeV where the phase 
transition occurs. The smaller role of the Coulomb field 
in the phase transition has recently been discussed exper- 
imentally in the framework of the Landau's description 
of phase transitions |17| . 

In conclusion, in this work we have addressed a gen- 
eral method for deriving densities and temperatures of 
fermions. The method has a validity O(j-) 3 and higher 



order terms might be included if needed. In the frame- 
work of the Constrained Molecular Dynamics model, 
which includes Fermi statistics, we have discussed colli- 
sions of heavy ions below lOOMeV/A and obtained densi- 
ties and temperatures at each bombarding energy. Know- 
ing the thermal energy of the system, we can derive the 
energy density and temperature reached during the col- 
lision. We have been able to bridge low energy phe- 
nomenology, i.e. particles evaporation from the surface, 
with the fragmentation of the system. Because of its 
general validity the approach could be applied to any 
fermionic system but for temperatures below the corre- 
sponding Fermi energies. We are also thinking about 
trapped Fermi gases [12| where the complete EOS could 
be derived following our method. Our approach is com- 
pletely at variance with previous ones based on classical 
mechanics. The results we have obtained here in a model 
case confirms that the classical approximation is unjus- 
tified. The tools we have proposed can be easily gener- 
alized to other fermion types, tritons, helions etc., and a 
comparative study of the EOS for different particles will 
be very interesting. We have seen in this work that differ- 
ent particles like neutrons and protons explore different 
density and temperature regions. Open problems such 
as Mott transitions, pairing etc. in low density matter 
might be addressed through a detailed study of the EOS 
using different fermions. A more conclusive study could 
be achieved if Boson-like particles could be included in 
the approach. This aspect will be the goal of our future 
work. 
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